clear


***REGRESSIONS

**TABLES 2-5
use "data_for_regression.dta", clear

drop if round==1
drop if region==1

macro def X fasting_mom log_duration log_interview_time interview_afternoon call_attempt_num  

gl E_FE "i.enumerator"
gl D_FE "i.interview_day"


**Table 2
xtreg dds_mom treat i.round, fe cluster(kebelle)
estimate store m1
estadd ysumm
xtreg dds_mom treat $X $Y i.round $D_FE $E_FE , fe cluster(kebelle)
estimate store m2
estadd ysumm
xtreg mindds_mom treat i.round , fe cluster(kebelle)
estimate store m3
estadd ysumm
xtreg mindds_mom treat $X i.round $D_FE $E_FE, fe cluster(kebelle)
estimate store m4
estadd ysumm
xtreg mindds_mom4 treat i.round , fe cluster(kebelle)
estimate store m5
estadd ysumm
xtreg mindds_mom4 treat $X i.round $D_FE $E_FE, fe cluster(kebelle)
estimate store m6
estadd ysumm
estout m1 m2 m3 m4 m5 m6, cells(b(star fmt(2)) ci(par fmt(2))) starlevels( * 0.10 ** 0.05 *** 0.010) stat(ymean r2 N, fmt(3 2 0))
esttab m1 m2 m3 m4 m5 m6 using "Table2.rtf",  keep(treat 3.round) b(%9.2f) ci(%9.2f) star(* 0.10 ** 0.05 *** 0.01) stat(ymean r2 N, fmt(3 2 0)) indicate("controls=$X" "interview day =*interview_day*" "enumerator=*enumerator*") coeflabels(3.round round _cons constant) title("Table 2: Effects on diet diversity of mothers") addnotes("Notes: controls included are a dummy variable indicating whether the mother was fasting, duration of interview, time of interview, a dummy variable if interview was conducted in the afternoon, and number of call attempts. Interview day and enumerator are day of the week and enumerator dummy variables, respectively") rtf replace nogaps
estimates drop _all 


**Table 3
xtreg staples_beans_nuts treat i.round, fe cluster(kebelle)
estimate store m1
estadd ysumm
xtreg staples_beans_nuts treat $X i.round $D_FE $E_FE, fe cluster(kebelle)
estimate store m2
estadd ysumm
xtreg animal_source treat i.round, fe cluster(kebelle)
estimate store m3
estadd ysumm
xtreg animal_source treat $X i.round $D_FE $E_FE, fe cluster(kebelle)
estimate store m4
estadd ysumm
xtreg veg_fruits treat i.round, fe cluster(kebelle)
estimate store m5
estadd ysumm
xtreg veg_fruits treat $X i.round $D_FE $E_FE, fe cluster(kebelle)
estimate store m6
estadd ysumm
estout m1 m2 m3 m4 m5 m6, cells(b(star fmt(3)) ci(par fmt(2))) starlevels( * 0.10 ** 0.05 *** 0.010) stat(ymean r2 N, fmt(3 2 0))
esttab m1 m2 m3 m4 m5 m6 using "Table3.rtf",  keep(treat 3.round) b(%9.2f) ci(%9.2f) star(* 0.10 ** 0.05 *** 0.01) stat(ymean r2 N, fmt(3 2 0)) indicate("controls=$X" "interview day =*interview_day*" "enumerator=*enumerator*") coeflabels(3.round round _cons constant) title("Table 4: Effects on diet diversity of mothers by food groups") addnotes("Notes: controls included are a dummy variable indicating whether the mother was fasting, duration of interview, time of interview, a dummy variable if interview was conducted in the afternoon, and number of call attempts. Interview day and enumerator are day of the week and enumerator dummy variables, respectively") rtf append nogaps
estimates drop _all 	


**Table 4	
xtreg dds_mom treat $X i.round $D_FE $E_FE if median_mother_age==0, fe cluster(kebelle)
estimate store m1
estadd ysumm
xtreg dds_mom treat $X i.round $D_FE $E_FE if median_mother_age==1, fe cluster(kebelle)
estimate store m2
estadd ysumm
xtreg dds_mom treat $X i.round $D_FE $E_FE if quint_education==1, fe cluster(kebelle)
estimate store m3
estadd ysumm
xtreg dds_mom treat $X i.round $D_FE $E_FE if quint_education==2, fe cluster(kebelle)
estimate store m4
estadd ysumm
xtreg dds_mom treat $X i.round $D_FE $E_FE if hhsize_hh<=5, fe cluster(kebelle)
estimate store m5
estadd ysumm
xtreg dds_mom treat $X i.round $D_FE $E_FE if hhsize_hh>5, fe cluster(kebelle)
estimate store m6
estadd ysumm
xtreg dds_mom treat $X i.round $D_FE $E_FE if quint_asset==1, fe cluster(kebelle)
estimate store m7
estadd ysumm
xtreg dds_mom treat $X i.round $D_FE $E_FE if quint_asset==2, fe cluster(kebelle)
estimate store m8
estadd ysumm
estout m1 m2 m3 m4 m5 m6 m7 m8, cells(b(star fmt(3)) ci(par fmt(2))) starlevels( * 0.10 ** 0.05 *** 0.010) stat(ymean r2 N, fmt(3 2 0))
esttab m1 m2 m3 m4 m5 m6 m7 m8 using "Table4.rtf",  keep(treat 3.round) b(%9.3f) ci(%9.3f) star(* 0.10 ** 0.05 *** 0.01) stat(ymean r2 N, fmt(3 2 0)) indicate("controls=$X" "interview day =*interview_day*" "enumerator=*enumerator*") coeflabels(3.round round _cons constant) title("Table 4: Heterogeneous effects on diet diversity of mothers by mother's age and education, household size, and baseline householdd wealth") addnotes("Notes: controls included are a dummy variable indicating whether the mother was fasting, duration of interview, time of interview, a dummy variable if interview was conducted in the afternoon, and number of call attempts. Interview day and enumerator are day of the week and enumerator dummy variables, respectively") rtf append nogaps
estimates drop _all 


**Table 5
xtreg dds_child treat i.round if age_child_months<36, fe cluster(kebelle)
estimate store m1
estadd ysumm
xtreg dds_child treat $X i.round $D_FE $E_FE if age_child_months<36, fe cluster(kebelle)
estimate store m2
estadd ysumm
xtreg mindds_child4 treat i.round if age_child_months<36, fe cluster(kebelle)
estimate store m3
estadd ysumm
xtreg mindds_child4 treat $X i.round $D_FE $E_FE if age_child_months<36, fe cluster(kebelle)
estimate store m4
estadd ysumm
estout m1 m2 m3 m4, cells(b(star fmt(3)) ci(par fmt(2))) starlevels( * 0.10 ** 0.05 *** 0.010) stat(ymean r2 N, fmt(3 2 0))
esttab m1 m2 m3 m4 using "Table3.rtf",  keep(treat 3.round) b(%9.2f) ci(%9.2f) star(* 0.10 ** 0.05 *** 0.01) stat(ymean r2 N, fmt(3 2 0)) indicate("controls=$X" "interview day =*interview_day*" "enumerator=*enumerator*") coeflabels(3.round round _cons constant) title("Table 5: Effects on diet diversity of children under 36 months") addnotes("Notes: controls included are a dummy variable indicating whether the mother was fasting, duration of interview, time of interview, a dummy variable if interview was conducted in the afternoon, and number of call attempts. Interview day and enumerator are day of the week and enumerator dummy variables, respectively") rtf append nogaps
estimates drop _all 

	
	
*TABLES S2-S4
use "data_for_regression.dta", clear

drop if round==2
drop if region==1

macro def X fasting_mom food_insecure  


**supplmentary table 2
xtreg dds_mom treat i.round, fe cluster(kebelle)
estimate store m1
estadd ysumm
xtreg dds_mom treat $X  i.round  , fe cluster(kebelle)
estimate store m2
estadd ysumm
xtreg mindds_mom4 treat i.round , fe cluster(kebelle)
estimate store m3
estadd ysumm
xtreg mindds_mom4 treat $X i.round , fe cluster(kebelle)
estimate store m4
estadd ysumm
estout m1 m2 m3 m4, cells(b(star fmt(2)) ci(par fmt(2))) starlevels( * 0.10 ** 0.05 *** 0.010) stat(ymean r2 N, fmt(3 2 0))
esttab m1 m2 m3 m4 using "Supplementary_table2.rtf",  keep(treat 3.round) b(%9.2f) ci(%9.2f) star(* 0.10 ** 0.05 *** 0.01) stat(ymean r2 N, fmt(3 2 0)) indicate("controls=$X") coeflabels(3.round round _cons constant) title("Supplementary Table 2: Effects on diet diversity of mothers") addnotes("Notes: controls included are a dummy variable indicating whether the mother was fasting and food security status of the househol") rtf replace nogaps
estimates drop _all 
	
	
**supplmentary table 3
xtreg staples_beans_nuts treat i.round, fe cluster(kebelle)
estimate store m1
estadd ysumm
xtreg staples_beans_nuts treat $X i.round , fe cluster(kebelle)
estimate store m2
estadd ysumm
xtreg animal_source treat i.round, fe cluster(kebelle)
estimate store m3
estadd ysumm
xtreg animal_source treat $X i.round , fe cluster(kebelle)
estimate store m4
estadd ysumm
xtreg veg_fruits treat i.round, fe cluster(kebelle)
estimate store m5
estadd ysumm
xtreg veg_fruits treat $X i.round , fe cluster(kebelle)
estimate store m6
estadd ysumm
estout m1 m2 m3 m4 m5 m6, cells(b(star fmt(2)) ci(par fmt(2))) starlevels( * 0.10 ** 0.05 *** 0.010) stat(ymean r2 N, fmt(3 2 0))
esttab m1 m2 m3 m4 m5 m6 using "Supplementary_table3.rtf",  keep(treat 3.round) b(%9.2f) ci(%9.2f) star(* 0.10 ** 0.05 *** 0.01) stat(ymean r2 N, fmt(3 2 0)) indicate("controls=$X") coeflabels(3.round round _cons constant) title("Supplementary Table 3: Effects on diet diversity of mothers by food groups") addnotes("Notes: controls included are a dummy variable indicating whether the mother was fastingnd food security status of the household") rtf append nogaps
estimates drop _all 	

	
**Supplementary table 4
xtreg dds_child treat i.round if age_child_months<36, fe cluster(kebelle)
estimate store m1
estadd ysumm
xtreg dds_child treat $X i.round  if age_child_months<36, fe cluster(kebelle)
estimate store m2
estadd ysumm
xtreg mindds_child4 treat i.round if age_child_months<36, fe cluster(kebelle)
estimate store m3
estadd ysumm
xtreg mindds_child4 treat $X i.round  if age_child_months<36, fe cluster(kebelle)
estimate store m4
estadd ysumm
estout m1 m2 m3 m4, cells(b(star fmt(2)) ci(par fmt(2))) starlevels( * 0.10 ** 0.05 *** 0.010) stat(ymean r2 N, fmt(3 2 0))
esttab m1 m2 m3 m4 using "Supplementary_table4.rtf",  keep(treat 3.round) b(%9.2f) ci(%9.2f) star(* 0.10 ** 0.05 *** 0.01) stat(ymean r2 N, fmt(3 2 0)) indicate("controls=$X") coeflabels(3.round round _cons constant) title("Supplementary Table 4: Effects on diet diversity of children under 36 months") addnotes("Notes: controls included are a dummy variable indicating whether the mother was fasting and food security status of the household") rtf append nogaps
estimates drop _all 
	


	